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Background and purpose — Previous studies of soft tissue sar- 
coma (STS) have identified a number of possible prognostic fac- 
tors; however, the majority of these include highly selected popu- 
lations, with unclear validation of data and insufficient statistical 
methods. We identified prognostic factors in a validated, popu- 
lation-based 30-year series of STS treated at a single institution, 
using an advanced statistical approach. 

Patients and methods — Between 1979 and 2008, 922 adult 
patients from western Denmark were treated at the Aarhus Sar- 
coma Center for non-metastatic STS in the extremities or trunk. 
The endpoints were local recurrence (LR) and disease-specific 
mortality (DSM). Prognostic factors were analyzed using a pro- 
portional hazard model, including continuous variables as cubic 
splines. Directed acyclic graphs were used to depict the causal 
structure. 

Results — The 5-year LR was 16% and the 5-year DSM was 
24%. Important prognostic factors for both LR and DSM were 
age, duration of symptoms, tumor size, grade, margin, and radio- 
therapy, while anatomical location (upper, lower extremity, trunk) 
was prognostic for DSM. 

Interpretation — In this population-based series of adult, 
non-metastatic STS, we included directed acyclic graphs, cubic 
splines, and a competing risk model in order to minimize bias, 
and demonstrated that these statistical methods are feasible. 
Using these statistical methods on a large, validated dataset, we 
excluded depth as a prognostic factor and established that age, 
duration of symptoms, size, grade, margin, and radiotherapy 
were important prognostic factors for both local recurrence and 
disease-specific mortality. 



Over the past 3 decades, numerous studies on soft tissue sar- 
coma (STS) have investigated prognostic factors for local 
recurrence and survival (Pisters et al. 1996, Kattan et al. 
2002, Stojadinovic et al. 2002, Zagars et al. 2003, Weitz et 
al. 2003, Gronchi et al. 2005, Gutierrez et al. 2007, Jebsen et 
al. 2008, Biau et al. 2011, 2012). While the prognostic value 
of some factors, e.g. grade and surgical margin, are generally 
accepted, the value of other factors such as age and location 
is still uncertain. 

To investigate this discrepancy between previous studies, 
we reviewed the literature and 3 core problems emerged: 
small sample sizes, heterogeneity of study populations, and 
differences in statistical methods used. 

The rarity of STS, with an incidence of 1.4 per 100,000 
(Maretty-Nielsen et al. 2013), makes a sufficient sample size 
challenging to obtain, and studies are often based on few 
patients. Thus, these studies may not have sufficient power to 
identify prognostic factors, or may identify them by chance. 

The heterogeneity of study populations in the literature — 
inclusion of selected subtypes only, high-grade tumors, or 
patients from centers with major tertiary referral practices and 
a larger proportion of large, high-grade, recurrent, or other- 
wise complicated sarcomas — results in possible selection bias 
and low generalizability. Few authors (e.g. Gustafson 1994) 
have used population-based data, and when they have, the 
validity of data sources has not been described or is unknown. 

The statistical approaches used to analyze data may intro- 
duce several problems, including a substantial overestimation 
of the failure risk if competing risk is not taken into account, 
biased estimates if the proper confounding variables are not 
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included in the adjusted analyses, and risk of over-optimized 
or irreproducible estimates and residual confounding if the 
data are analyzed using categorical variables instead of con- 
tinuous variables (Fine and Gray 1999, Royston et al. 2006, 
Shrier and Piatt 2008). Previously, statistical methods for han- 
dling these problems were not easy to use; however, recent 
implementations in statistical software have made it feasible 
to address these problems properly. 

The aim of this study was to identify prognostic factors in 
a validated, population-based series of adult STS, using an 
improved statistical approach with cubic spline regression in a 
competing risk setting. 



Patients and methods 

The Danish healthcare system provides tax-funded, free-of- 
charge healthcare for all residents. Since 1968, all the citizens 
of Denmark have been assigned a unique 10-digit civil per- 
sonal registration number (CPR number). These numbers are 
used in all Danish administrative registries and clinical data- 
bases. This allows linkage of databases at the level of the indi- 
vidual patient and ensures complete follow-up. 

Practice at the Aarhus Sarcoma Center (ASC) 

Since 1979, treatment of sarcoma in western Denmark (with 
2.5 million inhabitants) has been centralized at the Sarcoma 
Center of Aarhus University Hospital (ASC). The diagnostic 
program follows international guidelines and includes a clini- 
cal examination, diagnostic imaging, biopsy, and histopatho- 
logical evaluation conducted by an experienced multidisci- 
plinary team (Co-operative Cancer Departments 2007). In the 
study period, a 3 -tiered malignancy-grading system based on 
cellularity, necrosis, anaplasia, and number of mitoses was 
used (Jensen et al. 1991). 

The primary treatment is surgery aiming at a wide margin 
(including radical margins), according to the principles of 
Enneking et al.(1980). Low-grade tumors are mainly treated 
with marginal or wide margins (1 cm surrounding cuff of 
normal tissue), while deep intermediate and high-grade tumors 
are treated with wide margins and postoperative radiotherapy. 
The radiotherapy consists of fractions of 2 Gy to a total dose 
of 50-60 Gy. 

Patients are followed for a minimum of 5 years after the last 
treatment, with intervals ranging from 3 to 6 months. Follow- 
up visits include clinical examination and chest radiography 
(intermediate- and high-grade tumors only), with MRI scans 
for patients with deep-seated tumors in recent years. 

Study population 

Between January 1979 and December 2008, 1,753 consecu- 
tive patients were treated for STS at the ASC and registered 
in the Aarhus Sarcoma Registry (ASR) (Vraa et al. 1998). The 
ASR has recently been systematically validated by reviewing 
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Figure 1 . Flow chart of patients treated for a malignant soft tissue tumor 
at Aarhus Sarcoma Center in the period 1979-2008, with number of 
patients (n), reasons for exclusion, and the study population of adult 
patients with non-metastatic soft tissue sarcoma in the extremities or 
trunk wall. 



the medical files (gold standard) for all patients registered in 
the ASR using standardized forms. The completeness of the 
ASR was > 85% for western Denmark in 1979-2008, when 
compared with the Danish Cancer Registry at the patient level 
(Maretty-Nielsen et al. 2013). The ASR collects basic patient 
data, specific data on tumor characteristics and treatment, and 
data on follow-up, local recurrence, distant metastasis, and 
death. The completeness of registered data is > 99%. Patients 
were excluded if they had certain histological types, if they 
had tumors not located in the extremities or trunk wall, if 
they were younger than 15 years, or if they had metastases 
at diagnosis (Figure 1). The study population thus comprised 
922 adult patients with non-metastatic STS of the extremity 
or trunk. 

Prognostic factors 

Based on a literature search, possible prognostic factors were 
selected and prioritized before the data analysis. Based on 



Acta Orthopaedica 2014; 85 (3): 323-332 



325 



the number of events in the study (i.e. in order to have suf- 
ficient statistical power), the following factors were included: 
age at diagnosis, duration of symptoms, tumor size, location, 
depth, grade, surgical margin, radiotherapy, and calendar year 
of diagnosis. Directed acyclic graphs were used to depict a 
possible causal relationship between the prognostic factors 
selected, possible confounding variables, and the outcomes 
(Figure 2, see Supplementary data). 

Tumor size was recorded as the largest diameter (in cm), 
based on the fixed pathological specimen. However, if patients 
were not treated surgically or if the pathology report was 
insufficient, diagnostic images were used to determine size. 
Tumors located in the shoulder and gluteal area were clas- 
sified as trunk tumors. Margins were defined, based on the 
pathological specimen, as intralesional/marginal if the inci- 
sion was within the tumor or the pseudocapsule, or as wide if 
the tumor was surrounded by a cuff of normal tissue. 

Statistics 

Descriptive statistics are reported as median and proportions; 
they were compared using the Kruskal-Wallis test and the 
chi-squared test. The outcomes assessed were local recur- 
rence (LR) and disease-specific mortality (DSM). DSM was 
defined as death from sarcoma (including patients with locally 
advanced tumors) or death from metastatic disease. Patients 
were followed from the date of diagnosis (DSM analysis) or 
the date of their last primary treatment (LR analysis) to date 
of outcome, emigration, or end of the study period. The LR 
and DSM were estimated as cumulative incidence rates, with 
95% confidence intervals (CIs) based on the pseudo- value 
approach (Fine and Gray 1999, Graw et al. 2009). 

The possible prognostic factors were analyzed using a 
Cox proportional hazard model, estimating both crude and 
adjusted hazard ratios (HRs) with CI. The confounding vari- 
ables, as given in Tables 3 and 4 and Figure 2 (see Supplemen- 
tary data), were selected based on the approach described by 
Shrier and Piatt (2008). Age, duration of symptoms, and tumor 
size were analyzed continuously as 4-knotted cubic splines, 
with the twenty-fifth percentile as reference (Smith 1979). 
Adjusted cumulative incidence curves were constructed using 
the Stata command stcurv (Cleves 2000). The proportional 
hazard assumption was tested for all the variables included, 
both crude and adjusted, using log-minus-log plots and no 
contradictions were found. Competing risks were death from 
all causes in the analyses of LR and death from causes other 
than sarcoma in the analyses of DSM. Only patients who were 
macroscopically disease-free (i.e. with no macroscopic resid- 
ual tumor) after primary treatment (n = 879) were included in 
the analyses of LR. All tests were two-sided and any p- value 
< 0.05 was considered significant. Analyses were performed 
using Stata statistical software, version 13.1. 

Ethics 

The study was conducted in accordance with the Helsinki 



Declaration. It was approved by the Danish Data Protection 
Agency (2007-58-0010) and the Danish Health and Medicines 
Authority (7-604-04-2/262/KWH). 

Results 

Patient and tumor characteristics 

Overall, 922 adult patients in western Denmark were diagnosed 
with non-metastatic STS of the extremity or trunk wall from 
1979 through 2008. The median age was 60 (15-95) years and 
52% were males (Table 1). Over the study period, an increas- 
ing proportion of patients were referred without previous treat- 
ment and for primary diagnostics. At the same time, the pro- 
portion of patients treated with amputations decreased and the 
proportion of patients treated with radiotherapy increased. The 
proportion of wide resections decreased in the final 10-year 
period. The most frequent tumor types were malignant fibrous 
histiocytoma, liposarcoma, and leiomyosarcoma (Table 2). 
After primary treatment, 95% of the patients were macroscopi- 
cally disease-free. The median follow-up time was 7 (0-34) 
years. The median follow-up time for patients who were still 
alive at the end of the study period was 13 (3-34) years. 

Local recurrence (Figures 3-5 and Table 3) 

The 5 -year LR rate was 16% (CI: 13-18). The correlation 
between age and the risk of LR was approximately linear 
with a constant risk from 15 to 55 years, after which the risk 
increased successively. The relationship between the risk of 
LR and symptom duration was "J"-shaped, where an increase 
in the duration of symptoms in the first 12 months was associ- 
ated with a pronounced drop in the risk of LR, but after that 
the risk of LR increased steadily the longer the duration of 
symptoms. An increasing tumor size, for tumors < 5 cm or 
> 15 cm, was associated with an increasing risk of LR. In the 
confounder- adjusted analyses, age > 55 years (vs. 45 years), 
duration of symptoms < 3 months (vs. 3 months), tumor size 
5-7 cm or > 20 cm (vs. 4 cm), grade 3, and intralesional/mar- 
ginal excision were adverse prognostic factors, while duration 
of symptoms > 3 months and < 24 months (vs. 3 months), 
tumor size < 4 cm (vs. 4 cm), and treatment with radiotherapy 
were favorable prognostic factors. 

Disease-specific mortality (Figures 3, 4, and 6, and 
Table 4) 

The 5- and 10-year DSM rates were 24% (CI: 21-27) and 29% 
(CI: 26-32). The risk of DSM was nearly constant between 
15 and 45 years, after which it increased. Increasing duration 
of symptoms was associated with a decreasing risk within 9 
months, but after that it was approximately constant. Increas- 
ing tumor size was associated with an increasing risk, and this 
was most pronounced for tumors < 5 cm. Adverse prognostic 
factors were age > 45 years (vs. 45 years), tumor size > 4 
cm (vs. 4 cm), trunk and lower extremity location, grade 2-3 
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Table 1 . Characteristics of adult patients treated for a non-metastatic soft tissue sarcoma in the extremi- 
ties or trunk wall at the Sarcoma Center of Aarhus University Hospital according to calendar year of 
diagnosis (n = 922). Unless otherwise stated, values are number (%) 





-i 979-1 988 


1989-1998 


1 999-2008 


Total 


n-valup a 


Sex 












FpitiaIp 
i d i idic 


1 04 (49) 


1 47 (50) 


1 90 (46^ 

i \J\j yt\j j 


441 (48) 




Male 


110 (51) 


148 


223 (54) 


481 (52) 


0.6 


Acip dianno^i^ vpar*? 












Median (range) 


61 (15-95) 


57 (15-92) 


60 (15-93) 


60 (15-95) 


0.06 


Duration of symptoms, months 












Mpdian (rannp^ k 


9 (1-672) 


7 M-600 1 ! 


6 (1-732) 


6 (1-732) 


0.06 














None 


55 (26) 


1 35 (46) 


238 (58) 


428 (46) 




uiuijoy 


32 (15) 


34 M2^ 


32 (8^ 


98 (1 1) 




Surnprv in^uffiripnt c 

Ul UUI Vj II IOUI 1 lol^/l 11 


1 1 5 (54) 


120 (41) 


1 35 (33) 


370 (40^ 




Rurnprv QiiffiHpnt c 


12 (6) 


6 (2) 


8 (2) 


26 (3) 


< 0.001 


Rpa^on for rpfprral 

1 IwQOUl 1 l\Jl 1 CIUI lul 












Primary diagnostic 


1 84 (86) 


284 (96) 


397 (96) 


865 (94) 




Loral rprurrpnrp 

LUV/UI 1 wvul 1 ul IvC 


30 (14) 


11 (A\ 

i i y+) 


1 6 (A) 


57 (6) 


< 0.001 


Tumor size cm 












Median (range) ^ 


6 (1-30) 


6 (1-28) 


7 (1-35) 


7 (1-35) 


0.7 


Depth 












Superficial 


60 (28) 


72 (24} 


130 (31) 


262 (28) 




Deep 


1 54 (72) 


223 (76) 


283 (69^ 


660 (72^ 


0.1 


Grade 












1 


36 (17) 


44 (15) 


67(16) 


147 (16) 




2 


40 (19) 


55 (19) 


47 (11) 


142 (15) 




3 


129 (60) 


196 (66) 


292 (71) 


617(67) 




Unclassifiable 


9(4) 


0(0) 


7(2) 


16(2) 


0.001 


Surgical treatment 


207 (97) 


286 (97) 


401 (97) 


894 (97) 


1.0 


Type e 












Resection 


153 (77) 


226 (81) 


368 (92) 


747 (85) 




Amputation 


47 (24) 


54 (19) 


31 (8) 


132 (15) 


< 0.001 


Margin f 












Wide 


157 (79) 


223 (79) 


256 (64) 


636 (72) 


< 0.001 


Intralesional/marginal 


43 (22) 


59 (21) 


143 (36) 


245 (28) 




Radiotherapeutic treatment 


28(13) 


107 (36) 


161 (39) 


296 (32) 


< 0.001 


Chemotherapeutic treatment 


9(4) 


2(1) 


12(3) 


23 (2) 


0.03 


Macroscopically disease-free 


202 (94) 


283 (96) 


394 (95) 


879 (95) 


0.7 



a Based on the Kruskal-Wallis test and the chi-square test. 
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e Only patients treated with surgery (n = 894) included. 15 missing values. 
f Only patients treated with surgery (n = 894) included. 13 missing values. 



or unclassifiable, and intralesional/marginal excision, while 
duration of symptoms > 15 and < 30 months, tumor size < 4 
cm (vs. 4 cm), and treatment with radiotherapy were favorable 
prognostic factors. Depth was not a prognostic factor. 



Discussion 

In this population-based study, using systematically validated 
data, we found a 5-year LR rate of 16% and a 5-year DSM 
of 24%. Important prognostic factors for both LR and DSM 
were age, duration of symptoms, tumor size, grade, margin, 
and radiotherapy, while anatomical location was prognostic 
for DSM. 



Methodological considerations 

The main strengths of the present study include the large 
sample size, the population-based design with complete fol- 
low-up, and the statistical methods used to minimize bias. 

The validity of the ASR has been reported in detail in a pre- 
vious paper (Maretty-Nielsen et al. 2013). The overall com- 
pleteness of patient registration in the ASR is > 85%; how- 
ever, the completeness differed according to calendar year at 
diagnosis, residency, and age. The completeness was lower in 
patients diagnosed before 1994, in patients from the south- 
ern region of Denmark, and in patients over 60 years of age. 
These differences may introduce a selection bias if the asso- 
ciation between the factors investigated and the outcome is 
different for patients not registered in the ASR. Based on this, 
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Table 2. Distribution of non-metastatic soft 
tissue sarcomas in the extremities and trunk 
wall according to histological types in adult 
patients treated at the Aarhus Sarcoma Center 
from 1979 to 2008 (n = 922) 
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we would expect any selection bias to have 
been minor. 

The possible prognostic factors we inves- 
tigated were selected based on a literature 
search, and they were specified before the 
data analysis. When selecting the possible 
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Figure 4. Adjusted hazard ratios (HR (solid line)) with 95% confidence intervals (dashed 
lines) for local recurrence (panels A, C, and E) and disease-specific mortality (B, D, and F) 
according to age, duration of symptoms, and tumor size, based on Cox proportional hazard 
analyses. Adjustment covariates were selected based on Figure 2 (see Supplementary 
data); no covariates were included in the analysis of age; duration of symptoms was adjusted 
for age and grade; tumor size was adjusted for duration of symptoms and grade. 




Figure 3. Cumulative incidence rate (solid black line) of local recurrence (panel A) and disease-specific mortality (panel 
B) with 95% confidence intervals (dotted gray lines). 
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Figure 5. Adjusted cumulative incidence of local recurrence by grade (A), surgical margin (B), and radiotherapy (C) based on a Cox proportional 
hazard regression. Grade was adjusted for age and histological subtype; margin was adjusted for age, size, depth, location, compartmentalization, 
grade, and year of diagnosis; and radiotherapy was adjusted for age, depth, grade, margin, and year of diagnosis. 



Table 3. Crude and confounder-adjusted Cox proportional hazard analyses of possible prognos- 
tic factors for local recurrence in adult patients with non-metastatic soft tissue sarcoma in the 
extremity or trunk wall (n = 879) 
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a Confounding variables were selected based on the directed acyclic graph depicted in Figure 


2 (see Supplementary data): location was adjusted for histological type; depth was adjusted for 


duration of symptoms, size, and histological type; grade was adjusted for age and histological 


type; margin was adjusted for age, size, depth, location, compartmentalization, grade, and year 


of diagnosis; radiotherapy was adjusted for age, depth, grade, margin, and year of diagnosis; no 


adjustments were included in the analysis of year of diagnosis. 




b Only patients treated with surgery (n = 876) included. 3 missing values. 





confounding factors that should be included in the adjusted 
analyses, many studies have used variations of forward selec- 
tion, backward elimination, or simply inclusion of all factors; 
however, this can result in biased estimates with CIs that are 
too narrow and p-values that are too low. To our knowledge, 
this is the first study of prognostic factors in STS in which 



directed acyclic graphs have been 
used to depict causal relationships 
in order to select possible con- 
founding factors and thus minimize 
bias. 

The continuous variables were 
included as cubic splines (Smith 
1979). When evaluating the results 
from other studies, these variables 
are often analyzed categorically 
with one cut-off point or continu- 
ous linearly (Royston et al. 2006). 
As seen in Figure 4, neither of 
these approaches is adequate, with 
the categorical method resulting in 
residual confounding and the con- 
tinuous linear method being a false 
assumption. Furthermore, when 
selecting cut-points the method is 
often not specified, or the "optimal" 
cut-point method is used, leading to 
over-optimized and irreproducible 
estimates. Likewise, different cut- 
points are often used, impairing the 
generalizability. 

The WHO classification, which 
was used to classify STS in this 
study, has recently been updated, 
resulting in the change of some his- 
tological subtypes, e.g. malignant 
fibrous histiocytoma, now called 
undifferentiated pleomorphic sarcoma, and fibrosarcoma, 
which is becoming increasingly rare due to classification as 
other more specific types (Fletcher et al. 2002). The pathol- 
ogy of the STS patients reported in this study has not been 
reviewed, which may have caused misclassification when 
compared to the new WHO classification. However, a histo- 
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Figure 6. Adjusted cumulative incidence of disease-specific mortality by location (A), grade (B), surgical 
margin (C), and radiotherapy (D) based on Cox proportional hazard regression. Location was adjusted 
for histological subtype; grade was adjusted for age and histological subtype; margin was adjusted for 
age, size, depth, location, compartmentalization, grade, and year of diagnosis; and radiotherapy was 
adjusted for age, depth, grade, margin, and year of diagnosis. 



pathological evaluation according to the (at that time) relevant 
classification was performed at the time of diagnosis by an 
experienced sarcoma pathologist, for all patients referred to or 
diagnosed/treated at the Aarhus Sarcoma Center. Furthermore, 
during the study period only 2-3 pathologists performed these 
evaluations, ensuring that the classification was consistent 
throughout the study period. For this reason, misclassifica- 
tion would be expected to be minor and of no significance to 
our study, where several histological types were considered 
together. 

The comparability with other prognostic studies may have 
been complicated by the grading system used in the ASR 
(Jensen et al. 1991). This is primarily based on the number 
of mitoses, taking the cellularity, necrosis, and anaplasia into 
account, with some specific histopathological subtypes pre- 
specified, e.g. rhabdomyosarcoma and extraosseous osteosar- 
coma as grade 3. Currently, the most common grading sys- 
tems are the National Cancer Institute (NCI) system and the 
Federation Nationale des Centres de Lutte contre le Cancer 
(FNCLCC) system (Costa et al. 1984, Trojani et al. 1984). 
These grading systems are based on the same principles as in 



the system we used. STS patients 
at the ASC have been registered 
according to the FNCLCC system 
since January 2009, and based on 
the limited discrepancies in grad- 
ing that we have found before 
and after January 2009, we would 
expect comparability issues to be 
minor. 

We assessed the outcomes LR 
and DSM using a competing risk 
model, which was not done in 
most previous studies. In cases in 
which the competing event is fre- 
quent, not using a competing risk 
model can lead to overestimation 
of the failure outcome (Fine and 
Gray 1999). Establishment of the 
correct cause of death is essen- 
tial in DSM, and even though we 
used validated data on the cause 
of death, this can be problematic 
in elderly patients, where comor- 
bidity is frequent. 

During the study period, some 
changes in the treatment of 
patients at ASC were observed. 
Consistent with reports from 
other sarcoma centers, the pro- 
portion of patients treated with 
amputation decreased from 24% 
in the period 1979-1988 to 8% in 
the period 1999-2008 (Gronchi et 
al. 2005, Jebsen et al. 2008, Biau et al. 2011, 2012). Overall, 
it appears that there has been a shift in treatment regime, with 
a tendency towards "closer" surgical removal and adjuvant 
treatment with radiotherapy. However, this shift did not affect 
the LR rate or the DSM rate, as seen in the insignificance of 
calendar year of diagnosis in our analyses. 

Overall rates 

The LR and DSM rates presented in our study were in overall 
agreement with those in previously published studies (Gus- 
tafson 1994, Pisters et al. 1996, Stojadinovic et al. 2002, Zagars 
et al. 2003, Gutierrez et al. 2007, Jebsen et al. 2008, Biau et al. 
2011). Such studies have found 5 -year LR rates between 8% 
and 22%, whereas there tended to be higher LR rates previ- 
ously, which is possibly explained by the non-competing risk 
setting. The majority of previous studies found 5 -year DSM 
rates between 20% and 45%, as compared to our results of 
24%. The lower rates observed might also be explained by the 
fact that surgical treatment during the study period was per- 
formed by a small number of experienced oncology surgeons 
working in close collaboration. 



330 



Acta Orthopaedica 2014; 85 (3): 323-332 



Table 4. Crude and confounder-adjusted Cox proportional hazard analyses of possible prognostic factors for 
disease-specific mortality in adult patients with non-metastatic soft tissue sarcoma in the extremities or trunk 
(n = 922) 
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15.75 (5.78-42.93) 


Unclassifiable 
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50 
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20 
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296 (32) 
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1.28 (1.00-1.64) 


0.64 (0.48-0.86) 


Calendar year of diagnosis 
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214 (23) 


26 


33 
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295 (32) 


22 
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0.93 (0.68-1 .26) 
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1999-2008 


413 (45) 


24 


29 


0.93 (0.69-1 .25) 


0.93 (0.69-1 .25) 



DSM: disease-specific mortality; HR: hazard ratio. 

a Confounding variables were selected based on the directed acyclic graph depicted in Figure 2 (see Supple- 
mentary data): grade was adjusted for age and histological type; location was adjusted for histological type; 
depth was adjusted for duration of symptoms, size, and histological type; margin was adjusted for age, size, 
depth, location, compartmentalization, grade, and year of diagnosis; radiotherapy was adjusted for age, depth, 
grade, margin, and year of diagnosis; no adjustments were included in the analysis of year of diagnosis. 
b Only patients treated with surgery (n = 894) included. 13 missing values. 



Prognostic factors 

Age > 45-55 years was associated with higher LR and DSM. 
The important effect of age on DSM, when analyzed continu- 
ously, has not been described before. Two recent studies inves- 
tigated the prognostic value of age when analyzed using cubic 
splines: Gronchi et al. (2005) reported no independent effect 
on either LR or DSM, while Biau et al. (201 1) found a statisti- 
cally significantly higher LR rate with increasing age. Stud- 
ies have shown increased LR and mortality with increasing 
age, even when reporting DSM; however, the majority of these 
have reported age as a binomial categorical variable (Kattan 
et al. 2002, Weitz et al. 2003, Zagars et al. 2003, Gutierrez et 
al. 2007). 

Size, grade, and margin were important prognostic factors 
for both LR and DSM, which is supported by the previous 
literature (Pisters et al. 1996, Kattan et al. 2002, Stojadinovic 
et al. 2002, Zagars et al. 2003, Weitz et al. 2003, Gronchi et 
al. 2005, Jebsen et al. 2008, Biau et al. 2011, 2012). Unex- 
pectedly, we found a tendency of unclassifiable tumors being 
associated with higher mortality than even high-grade tumors. 
One explanation could be that this group represents the most 
undifferentiated sarcomas with high aggressiveness. 



Very few studies have analyzed duration of symptoms in a 
confounder-adjusted setting, and only as a categorical vari- 
able. Previous results have been highly contradictory; some 
studies found no prognostic effect of duration of symptoms 
(Rougraff et al. 2007), others found that short duration of 
symptoms was associated with increased mortality (Saithna 
et al. 2008, Urakawa et al. 2013), and finally others found 
that short duration of symptoms was associated with reduced 
mortality (Nakamura et al. 201 1). One possible explanation of 
these contractions might be that the correlation between dura- 
tion of symptoms and mortality is "J"-shaped, as an expres- 
sion of rapid increase in size of high-grade tumors — and thus 
short duration of symptoms — while low-grade tumors grow 
slowly and in the end may de-differentiate into higher grades. 
This is further supported by the findings of Maguire et al. 
(1994), namely that the effect of duration of symptoms on sur- 
vival was not linear. In accordance with this, we found that 
both curves were almost "J"-shaped, with higher LR and DSM 
rates at both ends of the symptom-duration spectrum. 

The literature regarding the prognostic value of depth on 
mortality is controversial. In contrast to the majority of pub- 
lished studies (Pisters et al. 1996, Kattan et al. 2002, Weitz et 
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al. 2003, Gronchi et al. 2005), our data did not support depth 
as a prognostic factor for DSM, which is supported by a few 
studies, including a recent study comparing the sixth and 
the seventh version of the American joint committee's stag- 
ing system, where depth is no longer included (Rydholm and 
Gustafson 2003, Maki et al. 2013). In line with our study, most 
studies have found that anatomical location is an important 
prognostic factor for mortality (Pisters et al. 1996, Kattan et 
al. 2002, Stojadinovic et al. 2002, Zagars et al. 2003, Gutier- 
rez et al. 2007). In concordance with the majority of published 
studies (Pisters et al. 1996, Weitz et al. 2003, Jebsen et al. 
2008), but not all (Stojadinovic et al. 2002, Biau et al. 2011, 
2012), we found that location and depth were not prognostic 
factors for LR. 

Most studies (Jebsen et al. 2008, Biau et al. 2011, 2012), 
like ours, found that radiotherapy reduces the LR rate; how- 
ever, the characteristics of tumors for which radiotherapy is 
beneficial vary. Kay tan et al. (2003) reported an effect only in 
high-grade tumors, while Jebsen et al. (2008) reported effects 
in low-grade tumors also. Very few studies have investigated 
the effect of radiotherapy on survival, and while some found a 
significant impact (Gronchi et al. 2005, Gutierrez et al. 2007), 
others found the opposite (Weitz et al. 2003). A retrospec- 
tive study involving 6,960 patients showed that radiotherapy 
improved survival for high-grade tumors, while no improve- 
ment was seen for low-grade tumors. (Koshy et al. 2010). 
Based on the local practice at the Aarhus Sarcoma Center, 
subgroup analyses were not obtainable in our study. 

Conclusion 

We conducted a large study including 922 consecutive adult 
non-metastatic STS patients, representing the "reality" seen 
in everyday clinics. Since results have a direct impact on 
future treatment strategies as well as the prognosis given to 
our patients, obtaining unbiased results is essential. In order to 
properly address a wide range of biases we included directed 
acyclic graphs, cubic splines and a competing risk model, and 
showed that these statistical methods are feasible. Using these 
improved statistical methods on a large, validated dataset, 
we excluded depth as an independent prognostic factor, and 
established that duration of symptoms and radiotherapy were 
important prognostic factors for DSM. Moreover, the method 
confirmed the importance of several prognostic factors for LR 
and DSM previously identified in studies with major tertiary 
referral practices. 

Supplementary data 

Figure 2 is available at Acta's website (www.actaorthop.org), 
identification number 6776. 
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